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Abstract 

In this paper we calculate the high-energy quark-quark scattering amplitude, first in the 
case of scalar QCD, using Fradkin's approach to derive the scalar quark propagator in 
an external gluon field and computing it in the eikonal approximation. (This approach 
was also recently used by Fabbrichesi, Pettorino, Veneziano and Vilkovisky to study the 
four-dimensional Planckian-energy scattering in gravity.) We then extend the results to 
the case of "real" (i.e. fermion) QCD, thus deriving again, in a rather direct way, the 
results previously found by Nachtmann. The abelian case (QED) is also discussed in the 
Appendix. 



1. Introduction 



High-energy scattering of strong-interacting particles attracted the interest of many 
physicists between the late 1950s and the early 1960s, when the first accelerators capable 
of producing protons in the GeV region were built. Before Quantum Chromo-Dynamics 
(QCD) was developed, essentially two relevant models (both originally based on potential 
scattering) were proposed to describe the physics of strong-interacting particles: the Regge 
pole model [||] and the droplet model of Yang and collaborators || . Although the physical 
pictures of these two models were surely interesting and had some partial experimental 
confirmations, they had not a strong theoretical foundation upon first principles. 

With the advent of quantum field theories, theoretical physicists became more and 
more involved in the study of the high energy behaviour of gauge field theories and in 
particular of Quantum Chromo-Dynamics, which is now generally believed to be the 
correct theory of hadrons and their interactions. A lot of work has been done within the 
framework of perturbation theory in order to find systematic procedures for extracting 
the large s, with t fixed, behaviour of each amplitude and for summing these contributions 
using a leading-log or eikonal approximation scheme |§ ||. 

Anyway there is a class of soft high-energy scattering processes, i.e. elastic scattering 
processes at high squared energies s in the center of mass and small squared transferred 
momentum t (that is s — > oo and t < s, let us say \t\ < 1 GeV 2 ), for which QCD 
perturbation theory cannot be safely applied, since t is too small. Elaborate procedures 
for summing perturbative contributions have been developed, even if the results are not 
able to explain the most relevant phenomena. 

A phenomenological analysis of this kind of soft hadronic reactions at high energies 
(in particular of the diflractive and elastic scattering of two hadrons) can be performed 
using the so-called Pomeron- exchanged model || . Even if it is nowadays accepted that 
the properties of the Pomeron come from the multi-gluon exchange among the various 
partons in the process, it is still an open problem how to quantitatively develop this 
picture in the framework of Quantum Chromo-Dynamics. Phenomenological analyses of 
the experimental results and a lot of theoretical investigations suggest that the Pomeron 
couples to the single partons within the hadrons, so that it is meaningful to study the 
parton-parton scattering amplitude in detail, and in particular the quark-quark and 
quark-antiquark amplitudes. 
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The first non-perturbative analysis, based on QCD, of these high-energy scattering 
processes was performed by Nachtmann in |J . He studied the s-dependence of the quark- 
quark (and quark-antiquark) scattering amplitude by analytical means, using a functional 
integral approach and an eikonal approximation to the solution of the Dirac equation in 
the presence of a non-abelian external gluon field. 

In this paper we propose a new approach to high-energy quark-quark (and quark- 
antiquark) scattering, based on a first-quantized path-integral description of quantum 
field theory developed by Fradkin in the early 1960s [0]. In this approach one obtains 
convenient expressions for the full and truncated-connected scalar propagators in an ex- 
ternal (gravitational, electromagnetic, ...) field and the eikonal approximation can be 
easily recovered in the relevant limit. Knowing the truncated-connected propagators, one 
can then extract, a la LSZ, the scattering matrix elements in the framework of a func- 
tional integral approach. We remind that this same method has been recently adopted in 
|| in order to study Planckian-energy gravitational scattering. We shall make directly 
use of some results contained in that paper, applying them to the case of interest to us, 
i.e. Quantum Chromo-Dynamics. 

This paper is organized as follows. In Sect. 2 we calculate the high-energy quark- 
quark scattering amplitude in the case of scalar QCD, using Fradkin's approach to derive 
the scalar quark propagator in an external gluon field and then computing it in the eikonal 
approximation. In Sect. 3 we extend the results to the case of "real" {i.e. fermion) QCD, 
thus deriving again, in a rather direct way, the results previously found by Nachtmann in 
||. The high-energy quark-quark scattering amplitude turns out to be described by the 
expectation value of two light-like Wilson lines, running along the classical trajectories 
of the two colliding particles. 

The abelian case is discussed in the Appendix, where we recover the well-known result 



for the eikonal amplitude of the high-energy scattering in QED [[J |T0[ |TT . 

A discussion on the validity of the various approximations used and on some calculation 
approaches proposed in the literature is included and concludes the paper. 
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2. The case of scalar QCD 



We begin, for simplicity, with the case of scalar QCD, i.e. the case of a spin-0 quark 
(described by the scalar field 0) coupled to a non-abelian gluon field = A^T a , T a 
(a = 1, . . . , N 2 — 1) being the generators of the Lie algebra of the colour group SU(N C ). 
We limit ourselves to the case of one single flavour. The Lagrangian is: 



4 ^ 



(2.1) 



where, as usual, D 11 = d^+igA^ is the covariant derivative. We then make use of Fradkin's 
approach to write the propagator for the scalar field as a functional integral of the first 
quantized theory [|l^]. This approach was also recently used by Fabbrichesi, Pettorino, 
Veneziano and Vilkovisky in |5J to study the four-dimensional Planckian-energy scattering 
in gravity. We shall use some results derived in ||, adapting them to the case of interest 
to us, i.e. high-energy scattering in scalar QCD. The extension to "real", fermion, QCD 
will be discussed in the next section. 

The Green function (Feynman propagator) for the scalar field <fi in the external metric 
g^u and in the external abelian field A^ admits a representation in terms of a functional 
integral over trajectories ||: 

G(y,x\g,A) = 

dv I [ y J-g(X)dX fl }5(x - X(0))6(y - X{y)) x 

(2.2) 



x exp 



[i/2) / dr{g tlv X i "X v + 2eA^X fl - m 2 ) 
Jo 



As explained in ref . ||13|| , one can employ the path-integral formulation of quantum me- 
chanics also when A^ is a matrix-valued potential of the form A^ = A^T a , where the 
matrices T a are the internal symmetry generators. As a consequence, the action S itself 
is a matrix with internal symmetry indices and, in order to have gauge covariance of the 
path-integral, it is necessary to define the exponential of the action by path-ordering. 
One can thus easily generalize the result (2.2) to the case of a non-abelian external gluon 
field N 1 : 

G(y,x\g,A) = 



dv 



_yJ-g(X)dX»]6(x - X(0))5(y - X(u)) x 



xPexp 



[i/2) / dT{ 9ilv X»X v ■ 1 
Jo 



2gA li X» 



m 2 ■ 1) 



(2.3) 



4 



where 1 is the N c x N c unity matrix in the colour space and "P" means " path-ordered 1 
along the "histories" X M (r). Moreover g(X) = det^^X)), so that g(X) = — 1 in a flat 
space-time. The propagator (2.3) is already connected, i.e. vacuum diagrams have been 
already divided out. Yet, in order to derive the scattering matrix elements following the 
LSZ approach, we need to know the on-shell truncated-connected Green functions: these 
are obtained from the connected Green functions by removing the external legs calculated 
on-shell. The truncated-connected propagator in the momentum space is defined as: 

G t (p,p'\A) = lim ^^-G(p,p'\A)^^- , (2.4) 

where G(p,p'\A) is the Fourier transform of G(x,y\A): 

G(p,p'\A) = J d A x J d 4 yexp[i(p'x - py)]G(x,y\A) , (2.5) 

with: G(x, y\A) =< T((fi(x)$(y) >a- In ref.|| it is shown how to evaluate the truncated- 
connected propagator (2.4) starting from eq.(2.3). It is also shown how to compute 
Gt(p,p'\A) in the so-called eikonal approximation, which is valid in the case of scattering 
particles with very high energy (E = p° ~ |p| ^> m) and small transferred momentum 
q = p' — p (i.e. \t\ <C E, where t = q 2 ). In this limit the functional integral can be 
evaluated by means of a saddle-point approximation in which the classical trajectory 
is computed to the lowest non-trivial order, that is a straight line (in the Minkowski 
space-time) described by a classical particle (of mass m) moving with 4-momentum p M : 

X » = W + V —v . (2.6) 

m 

One thus finds, after properly generalizing the results of || to the case of a non-abelian 
gluon external field, the following expression for the truncated-connected propagator in 
the eikonal approximation: 



G t (p,p'\A) -2E J d 3 be iqb (Vexp -ig J °° A„(b + pr)p^dr 



[2.7) 



As said before, q = p' — p is the transferred momentum. In "<i 3 6" one must not include 
the component of 6 M which is parallel to For example, if p^ ~ p'^ ~ (E, E, 0, 0), one 
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has that p + ~ p' + ~ IE and p_ ~ p'_ ~ 0, where the following general notation has been 
used for a given 4- vector 

A + = A + A 1 , A_ = A - A 1 . (2.8) 

In this case one must take: d 3 b = d 2 h t db-, where h t = (b y ,b z ) is the component of 
6 M in the transverse plane (y, z) (while the two light-cone coordinates, x + = t + x and 
x_ = t — x, are sometimes called "longitudinal -1 coordinates). 

The expression (2.7) for the truncated-connected propagator in the eikonal approxi- 
mation will be our starting point for deriving the scattering amplitude of two high-energy 
spin-0 quarks, in the limit s — > oo and t < s. We thus consider the elastic scattering 
process of two (scalar) quarks with initial 4-momenta p\ and P2 and final 4-momenta p[ 
and p' 2 : 

0(pi)0(p 2 ) - <P(p'i)<P(p 2 ) • (2-9) 

In the Center-of-Mass reference System (CMS), taking the initial trajectories of the 
two quarks along the x-axis, the 4-momenta pi, p 2 , p[ and p' 2 are given, in the first 
approximation, by: 

Pl ~p[~(E,E,O t ) , p 2 ~ p' 2 ~ (E,-E, t ) . (2.10) 

Using the LSZ reduction formula and a functional integral approach, one finds the follow- 
ing expression for the scattering matrix element (with the plane wave functions normalized 

as: <f> p (x) = exp(-ipx)): 

< 0(Pi)0(P2)l(^-l)l0(Pi)0(P2) >= 

= 4{< G t ( Pl ,p;|A)G t (p 2 ,p' 2 |A) > A + < (7 t (pi , p 2 1 A) (7 t (p 2 , Pi I ^4) >a} • (2.11) 

z <f> 

is the renormalization constant of the field <p: 4>r(x) = Z^ 1 ^ 2 <p(x). The expectation 
value < 0(A) >a of an arbitrary functional 0(A) of the gluon field A 11 is defined as: 



< 



0(A) > A = 1 1 [dA]0(A) exp [-i | rf 4 xFr^,] {det^Z^ + m 2 ]}' 1 , (2.12) 
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where Z is the partition function for scalar QCD: 

Z = J[dA][d(/)}{d^}exp i J d A xL(<p,<p\A) 



[dA] exp 



d 4 xFrF, 



afiu 



{det^Z^ + m 2 ]}- 1 



(2.13) 



The determinant in eqs.(2.12) and (2.13) comes from the integration over the scalar de- 
grees of freedom. In fact, the Lagrangian L(0, 4>\A) in (2.1) is bilinear in the scalar fields 
and 0^: 

J d 4 x [(7^0)^0 - m 2 0V] = - / d 4 x^ [d»D^ + m 2 ] . (2.14) 

Therefore the functional integral over the scalar fields and 0t is an ordinary Gaussian 
integral and can be performed in the standard way, to give (apart from an irrelevant 
constant) {det[ZW M + m 2 ]}" 1 . 

The first expectation value in eq.(2.11) corresponds to the t-channel scattering of 
the two particles, with squared transferred momentum equal to: (p± — p[) 2 = t <C s. 
The second expectation value corresponds instead to the -u-channel scattering of the two 
quarks. In other words, the squared transferred momentum, flowing from one quark to 
the other, is equal to: 

(pi-p' 2 ) 2 = u = Am 2 -s-t^-s, (2.15) 

in the limit we are considering: s — > oo with t, m 2 ^ s. In this high-energy limit the 
contribution coming from the second expectation value in eq.(2.11) is smaller by at least a 
factor of s, when compared with the first expectation value, and hence is negligible. One 
can be easily convinced of this by considering the Feynman diagrams of the process in 
the perturbation theory: the diagrams which correspond to the second piece in eq.(2.11) 
have intermediate gluons carrying a big squared transferred momentum u ~ — s, so that 
their propagators suppress the corresponding amplitude. Therefore, in our limit: 



< 0(P , 1 )0(P , 2 )|(S - l)|0( Pl )0(p 2 ) >~ -1 < GtiputilAfafatilA) > A . (2.16) 



The diffusion amplitude M/j =< f\M\i > is given, as usual, by: 

< 0^)0^)1(^-1)10(^)0^2) > = 

= i{27if5^\P fm - P m ) < 0(^)0(^)^10(^)0^) > 



(2.17) 
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where p n = pi + P2 is the total initial 4-momentum and P/j n = p[ + p 2 is the total final 
4-momentum. In the following we shall evaluate the expectation value: 



f(s,t) =< G t (p!, p[ \A)G t (p 2 ,p 2 \A) > A . 



(2.18) 



Let us indicate with Xi(t) and x 2 {t) the classical trajectories of the two colliding particles 
in the Minkowski space-time: 



xf(r) = 6f+pjfr , x^t) = b% + pt^T 



(2.19) 



By virtue of eq.(2.10) we can thus write the line-integral in the eikonal propagator for 
the particle 1 as: 

/+oo I r+oo "I r+oo 

A^xjdxt ~ -/ A^ Xl )dx l+ = - A_{x l+ M-Mt)dx 1+ . (2.20) 
-oo Z J — oo Z J — oo 

Similarly, for the particle 2: 

/+oo I /"+oo I r+oo 

A^(x 2 )dx^ ~ - / A + (x 2 )dx 2 - = - I A + (b 2+ ,x 2 -,b 2t )dx 2 - ■ (2.21) 
-OO ^ J — OO Z J —oo 



Therefore, if we define the two Wilson paths: 
^i(6 1 ) = W 1 (6 1+ ,6 1 _,b lt ) = Pexp 
W 2 (b 2 ) = W 2 (b 2+ ,b 2 _,b 2t ) = Pexp 



~7,9 / A_(xi + ,6i_,bi t )dxi+ 

2 J —oo 

/ A + (6 2 +,a:2_,b2t)rfx2_ 

Z J -co 



, (2.22) 



we obtain the following expression for the eikonal propagator (2.7) of the particle 1: 

d'budh.e^ / rf6 1+ — — b lt ) 
~4P / d%e mbl -— W 1 (b 1 ) . (2.23) 



We have used the fact that oo, &i_, bit) = 1; moreover qi-b\ + ~ 0, since pi_ ~ 

~ 0. We proceed in a similar way for the eikonal propagator of the particle 2, using 



8 



the fact that W 2 (b2+, — oo, h 2t ) = 1 and q 2+ b 2 - — (since p 2+ ~ p' 2+ ~ 0). We thus 
obtain that: 

G t (p 2 ,p' 2 \A) ~ 4E f d%e i ^-—W 2 (b 2 ) . (2.24) 
J ob 2 - 

We can substitute eqs.(2.23) and (2.24) into eq.(2.18), thus obtaining the following ex- 
pression for f(s, t): 

f(s,t) =< GtiputilAjGtfatilA) > A 

~ 1QE 2 { d\ { d%e^ + ^———— < W 1 (b 1 - b 2 )W 2 (0) > A . (2.25) 
J J ob 2 — 

We have made use of the translational invariance of the expectation value < W\(bi) 
W 2 (b 2 ) >a to write it as < Wiipi — 62)^2(0) >a- In order to isolate, inside f(s,t), 
the 4-dimensional <ie/£a-function coming from the conservation of the total 4-momentum 
(which we expect on the basis of eq.(2.17)), it will be helpful to define the new variables 
£ and z in the CMS: 

^ h + b 2 



2 

It is clear that: d% x d% 2 = d 4 T,d 4 z. Moreover: 



61 - 62 . (2.26) 



(91(9(9 (91(9(9 . . 

+ , ~oZ = o"ov : i ~a ' ( 2 - 27 ) 



db l+ 2.9E+ dz + ' db 2 ^ 2<9£_ dz_ 

Still making use of the translational invariance of the expectation value < W±(bi 
62)1^2(0) >a=< Wi(z)W 2 (0) >a, one can write, after some trivial algebra: 



f(s,t) = 1QE 2 J d 4 Ee i(<Zl+92)E J d A ze^- q2) " 



x 



x < Wi(z+,0,z t )W 2 (0,-z_,O t ) > A 

= (27r)V 4) (gi + q 2 ) ■ 8E 2 J d 2 z t J dz + J dz^e mz x 

x < m(z + ,0,z t )W2(0,-z_,O t ) > A . (2.28) 

Since gi_ ~ ~ ~ 0) and g 2 + — (p 2+ ~ p' 2 + — 0)> from the delta-i unction in 
front of eq.(2.28) one also derives that: q\ + = —q 2+ ~ and g 2 - = — qi- — 0. Therefore 
one can take the phase-factor exp(iqiz) out of the integrals in the light-cone variables z + 
and z_. Thus, defining: 

/(M) = ^(2vr) 4 5( 4 )(P /m - P in )/(s,t) , (2.29) 
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we finally find that (let us observe that s ~ 4E 2 and q\ + q 2 = (p[ + p' 2 ) — (pi + P2) = 

Pjin Pin) • 

f(s, t) = -lis J d 2 z t e^ Zt < [Wti+oo, 0, zt) - 1] x [W 2 (0, +00, t ) - 1] > A . (2.30) 

In the last equation we have put: 

q=Px—pf x = -qi ^ (0, 0, -q lt ) = (0, 0, q) . (2.31) 

That is: q = — q u - So that: qz ~ —q z t . The expression (2.30) gives the amplitude as an 
explicit function of the squared transferred momentum t = q 2 = — q 2 (therefore: t < !). 
More explicitly, putting in evidence the colour indices of the scalar quarks, the scattering 
amplitude can be written as: 

M f i =< M)^M)\M\^Mi(jP2) >= -^f{s,t) 



Z 2 

i 



,00 z 



2s J d 2 z t e 1 ^ < [WM) - iyw 2 (0) - l] kl > A ■ (2.32) 



Therefore, the high-energy (s —>■ 00 and i C s) scalar quark-quark scattering amplitude 
turns out to be the Fourier transform, with respect to the transverse coordinates z t , of 
the expectation value of two light-like Wilson lines, separated by z t = (0, 0, z t ): 



n+00 

WiC^sWiC+oM.Zt) = Pexp' 



W 2 (0) = W 2 (0,+oo,O t ) = Pexp 



/+00 
Apizt +pir)p1dr 
-00 

/+00 
A^(p 2 T)p^dr 
-00 



(2.33) 



The space-time configuration of these two Wilson lines is shown in [Fig.l]. This picture 
was originally proposed by Nachtmann in ref. (for the case of "real" fermion QCD, 
which we shall discuss in the next section). He derived the s-dependence of the quark- 
quark (or quark-antiquark) scattering amplitude using a functional integral approach 
and an eikonal approximation to the solution of the Dirac equation in the presence of a 
non-abelian external gluon field. 

From eq.(2.32) it seems that the s-dependence of the scattering amplitude is all con- 
tained in the kinematic factor 2s in front of the integral. Yet, as it was pointed out by 
H. Verlinde and E. Verlinde in $M, this is not true: in fact, one can easily be convinced 
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Fig. 1. The space-time configuration of the two light -like Wilson lines W\ 
and W2 entering in the expression (2.32) for the high-energy quark-quark 
elastic scattering amplitude. 



(for example by making a perturbative expansion) that it is a singular limit to take the 
Wilson lines in (2.32) exactly light-like. As suggested in [14|, one can regularize this sort 
of "infrared" divergence by letting each line to have a small time-like component, so that 
they coincide with the classical trajectories for quarks with finite mass m. In this way 
one obtains a logs-dependence of the amplitude, as expected from ordinary perturbation 
theory. We will return to comment on this at the end of the next section. 



3. The case of "real" (fermion) QCD 

We want now to extend the results obtained in the previous section to the case (more 
interesting from the physical point of view) of "real" fermion QCD: that is a non-abelian 
gluon field coupled to a spin-| quark. 

As before, we limit ourselves to the case of one single flavour. The QCD Lagrangian 

is: 

m^,A) = $(i-fDp - m)^ - \F%F a ^ , (3.1) 
where, as usual, = <9 M + igA^ is the covariant derivative. 

Proceeding exactly as in the previous section, we use the LSZ reduction formula and 
a functional integral approach to write the quark-quark scattering matrix element (with 



11 



the plane wave functions normalized as: ip p (x) = u p exp(—ipx), with u p u p = 2m) as: 

= ^2"{< ^a(pi)G , t (pi,pi|A)yW /9 (pi) ■ U 1 (p 2 )Gt(p2,p' 2 \A) k iUs(p 2 ) >A 

+ < u^ip'^GtipuP^^kiUpipi) ■ Ua(p'i)Gt(p2,p[\A)uu s (p2) >a} • (3.2) 

is the renormalization constant of the spinor field ip: i/jr(x) = Z^ 1 ^ 2 iIj(x). We have 
adopted the following relativistic normalization for Dirac spinors: 



u a {p)^up(p) = 2p tx 5 a p , u a (p)up(p) = 2m5 a(3 , 



(3.3) 



where a and (5 are spin indices. 

The expectation value < 0(A) > A of an arbitrary functional 0(A) of the gluon field 
A^ is now defined as: 

< 0(A) > A = i J [dA]0(A) exp [-1 1 rf 4 xFr^] det[i 7 ^ - m] , (3.4) 

where Z is the partition function for fermion QCD: 

Z = J[dA][di/>][di/>i]exp i J d 4 xL(4>,ijj\A) 



[dA] exp 



— - [ d A rF flu F 

J U XI a 1 afJiU 



det[i7 M D M - m] 



(3.5) 



The determinant in eqs.(3.4) and (3.5) comes from the integration over the fermion degrees 
of freedom. 

As explained in the previous section, in the high-energy limit we are considering, 
s — > oo with t, m 2 <C s, the contribution coming from the second expectation value in 
eq.(3.2) (w-channel scattering) is smaller by at least a factor of s, when compared with 
the first expectation value (t-channel scattering), and hence is negligible. Therefore, in 
our limit: 



< AM^kM) I (S- l)\^ j p(p l )lpis(p2) >^ 

- ^2 < «a(p'i)G' t (pi,pi|A) ij ^(p 1 ) • u 1 (p' 2 )G t (p2,p' 2 \A) k iu s (p2) >A ■ (3.6) 
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We must thus evaluate the expectation value: 



f(.S,t) =< Uaip'^GtipiiP'xlAjijUpip!) ■ U y (p' 2 )G t (p2,p' 2 \A) k iU 5 (p 2 ) >A ■ (3.7) 

The truncated-connected fermion propagator in the momentum space is defined as: 

G t (p,p'\A) = Km t—^G{p,p'\A) P -^ , (3.8) 

where G(p,p'\A) is the Fourier transform (see eq.(2.5)) of G(x,y\A) =< T(i[)(x)'i/;(y) >a- 
In eq.(3.8) we have used the notation: a = 7 M a M . As we shall see in the following, in the 
high-energy limit we are considering we can make the following replacement in eq.(3.7): 

Ua{p')G t (p,p\A)Uf3(p) > 

— > l-u a (p) Uf3 (p) ■ G { t scaL \p,p'\A) = 5 af3 ■ Gt aL) (p,p'\A) , (3.9) 

where we have indicated with Gj ''(p, p'\A) the truncated-connected propagator for a 
scalar [i.e. spin-0) quark in the external gluon field N 1 . Thus we approximately obtain: 

f(s,t) ~ 5 aP 8, s - < Gt cal \p l) p\\A) ll G^ al \p 2) p' 2 \A) M > A . (3.10) 

Let us see how this approximation is justified. Using the fact that u(p')(p' — m)u(p) = 
u(p')(p — m)u{p) = and the definition (3.8) of the truncated-connected fermion propa- 
gator, one can easily derive the following perturbative expansion: 



[p + q\Y — m A + is 



Uf,(p) + .. (3.11) 

(2) 



u a (p')G t (p,p'\A)up(p) = u a (p') [-ig^ A li (q l )} {l) up{p) + 
+ u a (p 

where we have used the following compact notation: 

[F(qi, ■■-,?»)](«) = 

Al 1 !Vqn i2-n)W\pi-p-q l -...-q n )F{q l ,... 1 q n ). (3.12) 



(2tt) 4 "J (2tt) 

Now, the key-point for the discussion (see, for example, ref. || and references therein) is 
that, in the high-energy limit we are considering, we have p ~ p' and fermions retain their 
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large longitudinal momenta during their scattering process. In other words, the exchanged 
vector mesons are allowed to carry only transverse momenta, so that the important phase- 
space region in the nth-order term of the perturbative expansion (3.11) has the property 
that qi and m are negligible when compared to p in the numerator (and similarly the terms 
of the form q^qj are negligible when compared to pqi in the denominator). Therefore, apart 
from a factor (?) n_1 coming from the n — 1 propagators, the numerator N(qi, . . . , q n ) of 
the nth-order term in (3.11) can be approximated as: 

N(q 1 ,...,q n ) = 

= u a (p')Y n (p + qi + ... + q n -i + m)^- 1 . . . ^ (p + q t + m^u^p) x 
x(-igA^ n (q n )) . . . (-igA^fa)) 

~ ^ a ^ u ^ ' (^"X 2 ^"" 1 ) • • • (2p w )(-^^n(?n)) • • • (-igAM) 

= 5 aP ■ (-^2^A Atn (g n ))(-^2^- 1 ^ n _ 1 (g n _ 1 )) . . . {-ig2p^ A^)) . (3.13) 

Apart from the (ie/ta-function in front, which simply reflects the fact that fermions retain 
their helicities during the scattering process in the high-energy limit, this is exactly the 
term we would have expected at the numerator of the nth-order term in the perturbative 
expansion of the truncated-connected scalar propagator G[ scaL \p, p'\A) in the high-energy 
limit (the denominators in (3.11) are already equal to the scalar case!). In fact, as reported 
in Table 1 (for scalar QED) and Table 2 (for scalar QCD), the factors (—ig2p fl ) in (3.13) 
just come from the quark-quark-gluon vertex of the scalar theory in the high- 

energy limit (when p ~ p'). 

Table 1: The quark-photon vertices in scalar QED. 



Vertex : 


Feynman rule : 


<t>i( P2 )4>( Pl )A» 


-ie(pi +p 2 T 




i2e 2 g^ 



One must observe (see again Table 1 and Table 2) that the scalar theory has also an- 
other type of vertex of the form ^(pA^A 1 " (quark-quark-gluon-gluon). Yet one can easily 
be convinced that the contributions, due to these additional scalar vertices, to the nth- 
order term in the perturbative expansion of the truncated-connected scalar propagator 
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Table 2: The quark-gluon vertices in scalar QCD. 



Vertex : 


Feynman rule : 




-ig{Pl+P2Y(Ta)ij 




ig 2 g^{T a ,T b } i:i 



Gt(p,P'\A) in the high -energy limit are suppressed with respect to the contribution 
corresponding to (3.13), i.e. that coming only from the quark-quark-gluon couplings. 
This is essentially due to the fact that the quark-quark-gluon-gluon vertex does not 
carry momentum (see Table 1 and Table 2); moreover, when a vertex of this kind (of 
order 0(g 2 )) is inserted in place of two vertices of the quark-quark-gluon kind (each of 
them of order 0(g)), a free fermion propagator (almost on mass-shell) desappears. 

One can also verify directly that the quark-quark-gluon-gluon vertices do not bring 
contributions to the quark-quark elastic scattering amplitude at the leading order in the 
high-energy limit. For example one can evaluate the scattering amplitude up to the order 
0(g 4 ) in the scalar theory. The diagrams with only quark-quark-gluon vertices give rise 
to the following contribution (see, for example, ref.0 and references therein): 

M (S 4) ~ <Al - a(t) Ins] ■ {Gi) ijM + ig A sI(t) ■ (G 2 ) ijM , (3.14) 

v 

where G\ and G 2 are colour factors for the group SU(N C ): 

(Gi)ij t ki = (T a )ij(T a ) k i , (G 2 )ij t ki = (T a T b )ij(T a T b ) k i , (3.15) 
and the functions a(t) and I(t) are defined as: 

a{t) = -^ m , , (()= /|^_L_ (t= _ A .). (3,6) 

Instead one finds that the contributions coming from those diagrams with at least one 
quark-quark-gluon-gluon vertex are at most of order 0(s°) = 0(1) and so they are 
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suppressed with respect to the terms in (3.14). Therefore eq.(3.14) is just the quark- 
quark elastic scattering amplitude at the leading order in the high-energy limit, up to the 
perturbative order 0(g 4 ) in the scalar theory. 

Therefore we have found that, in the high-energy limit we are considering, the re- 
placement (3.9) is fully justified and we are left with the expression (3.10) for f(s,t). In 
the expectation value (3.10) we must consider G[ scaL \p, p'\A) as a functional of the gauge 
field equal to the truncated-connected propagator for a scalar (i.e. spin-0) quark in 
the external gluon field A^. At this point we can use the results of the previous section, 
where Gf Scal '\p, p'\A) was evaluated in the eikonal approximation (see eq.(2.7)). We then 
proceed exactly as in the previous section, so that we finally find the following expression 
for the high-energy quark-quark elastic scattering amplitude in (fermion) QCD ||: 



■ S aP S lS ■ 2s J d 2 z t e l(i Zt < [W^zt) - iyW 2 (0) - 1] H > A . (3.17) 



The notation is the same as in the previous section. 

In a perfectly analogous way one can also derive the high-energy scattering amplitude 
in the case of the abelian group U(l) (QED). The resulting amplitude is equal to eq.(3.17), 
with the only obvious difference that now the light-like Wilson lines W\ and Wi are 
functionals of the abelian field A^ (so they are not matrices). Thanks to the simple 
form of the abelian theory (in particular to the absence of self-interactions among the 
vector fields), it turns out that it is possible to explicitly evaluate (at least in the quenched 
approximation) the expectation value of the two Wilson lines: the details of the calculation 
are reported in the Appendix and one finally recovers the well-known result for the eikonal 
amplitude of the high-energy scattering in QED || |TU| flTTfl . 

We want to conclude by making a few comments on eq.(3.17) for the case of QCD. As it 
was pointed out at the end of the previous section, it is a singular limit to take the Wilson 
lines in (3.17) exactly on the light-cone. It is expected [ITJ] that a proper regularization 
of these singularities will give rise to the logs-dependence of the amplitude, as obtained 
by ordinary perturbation theory || jp and as confirmed by the experiments on hadron- 
hadron scattering processes. 

Therefore, the direct evaluation of the expectation value (3.17) is a highly non-trivial 
matter, being also strictly connected with the ultra-violet properties of Wilson-line oper- 
ators. This last subject was considered for the first time in ref. [O many years ago. More 
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recently, in ref . JIB) , it has been found that there is a correspondence between high-energy 
asymptotics in QCD and renormalization properties of the so-called cross singularities of 
Wilson lines. The asymptotic behaviour of the quark-quark scattering amplitude turns 
out to be controlled by a 2 x 2 matrix of the cross anomalous dimensions of Wilson lines. 

An alternative non-perturbative approach for the calculation of the expectation value 
(3.17) has been proposed in ref.|[L7||. It consists in studying the Regge regime of large 
energies and fixed momentum transfers as a special regime of lattice gauge theory on an 
asymmetric lattice, with a spacing do in the longitudinal direction and a spacing a t in the 
transverse direction. In the limit a^/at — > one gets a possibility to study correlation 
functions such as (3.17) and to obtain the desired logs-dependence of the amplitude. 

It would be interesting to investigate about the possibility to use asymmetric lattice 
gauge theory to measure directly the amplitude (3.17) by Monte Carlo simulations. At 
the moment, the only non-perturbative numerical estimate of (3.17), which can be found 
in the literature, is that of ref. (where it has been generalized to the case of hadron- 
hadron scattering): it has been obtained in the framework of the model of the stochastic 
vacuum. It is still an open question whether a more conventional approach based on 
numerical simulations of lattice QCD can or cannot be of some help. Some progress in 
this direction is expected in the near future. 
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Appendix: The abelian case. 



In this appendix we shall discuss the abelian case (See also refs. |TB| and [fT7|). The 
fermion-fermion (or scalar-scalar) electro-magnetic scattering amplitude, in the high- 
energy limit s — > oo and i C s, can be derived following exactly the same procedure 
used in Sect. 2 and Sect. 3 for the non-abelian case. The resulting amplitude (for the 
fermion case) is formally identical to eq.(3.17), with the only obvious difference that now 
the Wilson light-like lines W\ and W 2 are functions of the abelian field (so they are 
not matrices): 



— --7PT- <W<5 7 <5 ' 

:— >oo X , 



2s J d 2 z t e tcl Zt < [Wxizt) - l][W 2 (0) -1}> A . (A.l) 
The electro-magnetic light-like Wilson lines W\ and W 2 are defined as in (2.33): 

/+00 
-00 
r+00 

W 2 (0) = W 2 (0,+oo,O t ) = exp -ie / A^r^dr , (A.2) 

J —CO 

where z t = (0, 0, z t ) and e is the electric coupling-constant (electric charge). Thanks to the 
simple form of the abelian theory (in particular to the absence of self-interactions among 
the vector fields), it turns out that it is possible to explicitly evaluate the expectation 
value of the two Wilson lines, thus finally recovering the well-known result for the eikonal 
amplitude of the high-energy scattering in QED (see refs. [f|, [[U| and fLlp. Let us see 
the details of the calculation. 

We shall evaluate the expectation value < Wi(z t )W2(0) >a in the so-called quenched 
approximation, where vacuum polarization effects, arising from the presence of loops 
of dynamical fermions, are neglected. This amounts to setting det(iT[A]) = 1, where 
K[A] = ij^D^ — m is the fermion matrix. We shall also put ~ 1. Thus we can write 
that: 

< Wi(z t )W 2 (0) > A ~ i JidA^W^z^iO) , (A.3) 

where Sa = — \ J d 4 xF llu F^ v is the action of the electro-magnetic field and Z = J[dA fl ]e' lSA 
is the pure-gauge partition function. We then add to the pure-gauge Lagrangian L A = 
— ^F flu F flu a gauge-fixing term Lqf = — 2q(^ m ^m) 2 (covariant or Lorentz gauge). The 
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expectation value (A. 3) becomes, denoting Lq = La + Lqf'- 
<W 1 (z t )W 2 (0) >A ~-^J[dA]x 



x exp 



/ r+oo r+oo \ 

d xL Q -e A^(z t + p^pldT - e / A^p 2 r)p%dT 
J— oo J —oo / 



, (A.4) 



with Z' = J[dA] exp(i / d^xL^). It will be helpful, for the computation, to write the sum 
of the two line integrals as a 4-volume integral of the form: 

/+oo r+oo r 

A^zt + p x r)f x dr - e / A^p 2 r)p^dr = / d^xA^x) J"(x) , (A.5) 
-oo J —CO J 



where J tl (x) is a 4-vector source defined as: 



J^x) = Jf +) (x) + Jl^x) = -e[ef +) 5(x_)5(x t - z t ) + ef.^+W**)] , 



(-r 



: (+) 



(A.6) 



with: e^ + ) = (1,1, 0, 0) and = (1, —1, 0, 0). Therefore the exponent in eq.(A.4) can be 



written as: 



% J d 4 x{L% + A^) 



(A.7) 



After integration by parts, the integral of L$ can be written as an integral of a quadratic 
form of the electro-magnetic field A^\ 



J d*xL% = J d A x 



-\F^-^A,f 



I 



d A x 



--A,K^A U 



(A.8) 



and the (invertible!) matrix is so defined: 



= -g, v d 2 + (l - £) d,d u . 

The functional integral 

Z [J] = J [dA\ exp i J d A x(L% + A^r) 

= J [dA\ exp i J d 4 x [- 1 -A ii K iiV A v + A^ 



(A.9) 



(A.10) 
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can be evaluated with standard methods (completing the quadratic form in the exponent). 
One thus obtains: 



Z [J] = Z> • exp [1 J d A x J d A yJ»{x)D, u {x - y)J"(y) 



(A.ll) 



where Z' = Z [J = 0] = J[dA] exp(i J d A xL^) and is the inverse of the operator K^: 



J d A zK^[x - z)g xp D pv {z -y) = g^(x - y) 



(A.12) 



In other words, Dn V is the free photon propagator and is equal to: 



D^(x 



v) = J 



d 4 k e-^-v) ( k^K' 
9nv - {I - a) 



(2vr) 4 k 2 + ie 



k 2 



(A.13) 



In the following we shall choose the gauge-fixing parameter a equal to 1. From eqs.(A.4), 
(A.ll) and (A.13), we derive the following expression for the expectation value of the two 
Wilson lines (A.2): 



< Wi(^)W 2 (0) > A - ^7 / [dA]e i $ dix{L » +A ^ ) 

Zj J 



= Z [J]/Z' = exp 



d A k e - ik( - x - y) 



(2tt) 4 k 2 + ie 



(A.14) 



Proceeding in exactly the same way, we can also derive (always in the quenched approxi- 
mation) the expectation values of the single Wilson lines Wi(z t ) and W 2 (0): 



< Wi{z t ) > A ~ exp 

< W 2 (0) > A ~ exp 



d 4 k e~ ik ^-y) 



J d*x J d*yJ» +) {x)J {+) »{y) j (2?r)4 k2+ . £ 
ft r d 4 k e - ik ( x -y) 

J dAx J d4 y J (-)(*) J (-M J (27r)4 k2 + i£ 



, (A.15) 



where the 4-vector sources and are defined in eq.(A.6). Using the fact that: 
e (+) e (+)M = e ( l _) e (-)M — 0, one finds that: 



< W 1 (z t ) > A =< W 2 (0) > A ~ 1 



(A.16) 
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Finally, using the explicit form (A. 6) of the 4-vector source J^(x) to evaluate the double 
integral in eq.(A.14), one finds the well-known eikonal form of the high-energy electro- 
magnetic scattering amplitude J| [[K]] [p"Tf): 

M fi ~ -i • 5 af3 5 yS ■ 2s J d 2 z t e^ < [W x {zt) - 1] ■ [W 2 (0) - 1] > A 

~ -i ■ 5 aP 8 lS ■ 2s J d 2 z t e^{< W 1 (z t )W 2 (0) > A -1] 

~ -i ■ S a ^ s ■ 2s J d 2 z t e l * zt (exp -ie 2 j * 2 ^ _ - - lj . (A.17) 

We could have also proceeded in a slightly different way, by observing that the free 
photon propagator has a rather simple expression in the coordinate space. Explicitly, in 
the Lorentz gauge with a = 1: 



9nu 



Att 2 (x 2 — ie) 



(A.18) 



After inserting this expression into eq.(A.ll) and integrating with respect to x and y, we 
are left with the following expression for the scattering amplitude: 



M fi ~ -i ■ 5 al3 5 lS ■ 2s / d z z t e l * zt exp 



dad/3 



a(3 + z\ + ie 



(A.19) 



This expression is, of course, perfectly equivalent to that in eq.(A.17), as one can eas- 
ily verify by explicitly evaluating the two integrals in the exponents. In fact, the two 
exponentials in eq.(A.17) and eq.(A.19) turn out to be equal to: 



e ■ exp ( i— In |z| 



(A.20) 



where A is an infinite constant phase and is therefore physically unobservable. The origin 
of this infinite constant phase resides in the well-known fact that the fermion-fermion 
scattering amplitude in QED has infrared divergences, due to the emission of low-energy 
massless vector mesons. The traditional way to handle these infrared divergences is to 
introduce an infrared cutoff in the form of a vector meson mass fi. In this way the integral 
over k t in the exponent of eq.(A.17) is substituted by the following expression: 



d 2 k t e ~ ikt Zt 1 



(2tt) 2 kf + /i 2 2tt 



7T K oO*\*t\) 



(A.21) 
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where K is the modified Bessel function. In the limit of small \x this last expression can 
be replaced by: 

^» ( " |z|) ,2.-^ ln (r >|z ' 1 ) ' < A ' 22 > 

Absorbing |e 7 in \x and putting A = —(l/2n) ln(|e 7 ii) (so that A — > oo when /x — * 0), we 
just obtain the expression (A. 20) for the exponential in eq.(A.17). The integral in (A. 17) 
with respect to the transverse coordinate z t can be explicitly performed and one finally 
gets: 

where: a = e 2 /4n. This is just the standard eikonal formula for the fermion-fermion 
electro-magnetic scattering amplitude. 
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